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It is known from the solution of the two-body problem that an anisotropic dipolar interaction can 
give rise to s-wave scattering resonances, which are named as dipolar interaction induced resonaces 
(DIIR). In this letter, we study zero-temperature many-body physics of a two-component Fermi 
gas across a DIIR. In the low-density regime, it is very striking that the resulting pairing order 
parameter is a nearly isotropic singlet pairing and the physics can be well described by an s-wave 
resonant interaction potential with finite range corrections, despite of the anisotropic nature of 
dipolar interaction. The pairing energy is as strong as a unitary Fermi gas nearby a magnetic 
Feshbach resonance. In the high density regime, the anisotropic effect plays an important role. We 
find phase transitions from singlet pairing to a state with mixed singlet and triplet pairing, and then 
from mixed pairing to pure triplet pairing. The state with mixed pairing spontaneously breaks the 
time-reversal symmetry. 



With the development of STIRAP technique [1 , it is 
now very promising to achieve degenerate gases of polar 
molecules within few years. The dipole moment of a po- 
lar molecule can be tuned continuously by increasing the 
strength of applied external electric field. Because of the 
large permanent dipole moment, the dipolar interaction 
energy can be on the same order of magnitude as the 
Fermi energy. Previously, many studies of dipolar Fermi 
gases focus on p-wave pairing due to the anisotropic na- 
ture of dipolar interactions [2H1] ■ However, the transition 
temperature for p-wave superfluid is usually suppressed 
by the centrifugal barrier. Thus, if the strong dipolar in- 
teraction can also cause significant effects in the s-wave 
channel, it must be first seen in experiments once a chem- 
ically stable degenerate polar molecular gas is realized. 

The interaction potential between two molecules with 
dipole moment dz is given by (we also set H = fcs = 1 
throughout this letter): 



V d (r) = 



2D 1-2, cos 2 0, 



(1) 



where D = md 2 /2 is dipole length, r is the polar an- 
gle between r and z. At first, one can easily see that 
(s|Vd(r)|s) = for any s-wave wave- function |s). Thus, 
within Born approximation the interaction effect in the 
s-wave channel is always weak. However, because Vd(r) 
couples s-wave |s) to d-wave \d) 1 and through a second- 
order process, it induces an effective attractive potential 
proportional to —D 2 /r 4 in the s-wave channel [5]. Thus, 
as D increases, the potential becomes deeper and deeper, 
and therefore s-wave bound states will appear sequently. 
This will lead to a series of scattering resonances in the 
s-wave channel [51-110) . which are named as dipolar in- 
teraction induced resonances (DIIR) . Two-body solution 
shows that nearby an s-wave DIIR, the scattering am- 
plitude is dominated by s-wave component and is very 
isotropic for low-energy scattering [8] , and the scattering 
phase shift is strongly momentum dependent, which gives 
rise to a large and positive effective range proportional 



to D [ID] . This is strongly in contrast to a wide Feshbach 
resonance studied extensively before, where the effective 
range is negligible. 

One major difficulty in studying many-body physics 
of dipolar gases is that the bare interaction Eq. ([T]) has 
short range divergency. This is because the Fourier trans- 
form Vd(k) of Vd(r) approaches a constant as k — > oo, 
thus, it has the similar ultraviolet divergency as a S func- 
tion contact potential. One common approach to elimi- 
nate this divergency is to relate Vd(k — k') to scattering 
amplitude r^k' via Lippman-Schwinger equation as 



k.k' 



I/ d (k-k')-5> k , q ™y d (q-k'). (2) 



This approach is widely used in regularizing a ^-function 
interaction potential [11]. Nevertheless, for anisotropic 
dipolar interaction potential, such an equation is diffi- 
cult to solve. So far most of previous studies of many- 
body physics of dipolar Fermi gas are based on the Born 
approximation, which assumes Tk.k' = Vd(k— k'). It 
greatly simplifies the solution, but the resulting Tk.k' is 
always anisotropic and clearly it fails near a DIIR. On 
the other hand, so far the studies of DIIR are limited to 
either two-body physics [SMTP] or the high-temperature 
regime [TU] , 

In this letter we study the zero temperature many- 
body physics for a two-component Fermi gas across a 
DIIR. Here we choose an alternative way to overcome the 
regularization problem, that is, by introducing a regular 
model potential at short distance, which is given by 

Tr , , 2D 3 cos 2 T -1„, , , 

V(r) = 3* F(r/r a ) (3) 

m r 3 

where F(x) is chosen as F(x) — e~ x (x 3 /6 + x 2 /2 + x + 
1) — 1. We choose this particular form of model poten- 
tial because its Fourier transform V(k) has an analytical 
expression 



V(k) = 



8ttD3cos 2 k - 1 
IrrT (1 + fc 2 r 2 ) 2 



(4) 
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As shown in Fig. 1, F(x) — > —1 exponentially as x 3> 1, 
thus V(r) behaves the same as the bare dipole interac- 
tion Vd(r) when r 3> r Q , and F(x) —> X clS X — > 0, which 
ensures V(k) decays sufficiently fast as k — > oo. More im- 
portantly, we should note that only the physics nearby a 
DIIR is independent of how we implement the regular- 
ization. Thus, we need to first locate the DIIR for this 
particular potential by solving the two-body Scrodingcr 
equation. Indeed, we find the first DIIR at D/r = 8.81 
with a positive effective range r e = 1.15D, as shown in 
Fig. 1(b) and (d). Our following many-body study will 
be restricted to the regime nearby this resonance. 

Mean-field Theory. The many-body Hamiltonian for 
this two-component Fermi gas is given by 

H = E^LV'ak + w A U (W(k-k')n a/3 (k',q)(5) 

ak kk'q 

where £k = k 2 /(2m), a and j3 denote spin f and i, 
and n a/3 (k, q) = -0 Qq /2-kV' / 3q/2+k- Here the spin refers 
to hyperfine spin degree of a polar molecule |12) . while 
the dipolar interaction is originated from the electronic 
dipole moment, thus we can take the dipole interaction as 
spin independent. Therefore, this Hamiltonian possesses 
an SU (2) spin rotational symmetry. Following the stan- 
dard mean-field approach, we decouple the interaction 
term in the Cooper channel and obtain 

H ml = ^e k £ k 4 k +^^[nt /3 (k,0)A^(k)+h.c] 

ak a /3k 

- \ J2 A* /J (k)V- l (k,k')A o;9 (k') (6) 

a^k,k' 

where A Q(3 (k) = £ k , V(k-k')(n a/3 ( k ', 0)) and 
V-^V) satisfies £ k „ V~ l (k, k")V(k" - k') = <5 k>k ,. 
The four superfluid order parameters A a p can be 
organized into one singlet and three triplet compo- 
nents as A s = i(A n - A ;t ), AJ = ^(A n + A it ), 

Af = -^(A tt - A u ) and A? = ~^(A n + A u ). 
Tracing out the fermion part, it yields the free energy: 

Q °=l E [ek-2Tln(2cosh^)] 

k,j=l,2 

-^E[ A ^ 1 (k,k')A sk ^+ Yl KZVf 1 ^*!)^ 

k.k' u=x,y,z 

(7) 

where £ k = e k - /x, V a , t (k,k') - [V(k + k') ± 
F(k-k')]/2, £ jk = y^k — A| k . and A 2 2k = [L x ± 
V^Lf - L 2 ]/2 with Lx = E lA^I 2 + |A sk | 2 and 

L 2 = | E ( A "k) 2 " ( A S k) 2 |- Both Li and L 2 are 

v— x,y,z 

invariant under an SU (2) spin rotation. Therefore, in 




FIG. 1: (a) Regularized dipolar interaction potential V(r), 
and the dashed line indicates the long range dipolar part, (b- 
d) s-wave scattering length a s (b) p-wave scattering length a p 
(c) and the s-wave effective range r c (d) around a DIIR. The 
plot is around the first DIIR located at D/ro — 8.81. 

analogy with spin-1 Bose condensate [13], we can classify 
the triplet components into two different types similar 
as polar phase (with A| k = A t k, A^ k = A^ k = 0) and 
ferromagnetic phase (with A^ k = 0, A^ k = zA^ k = A t k). 
It is easy to show that the ferromagnetic configuration is 
always energetically unfavored and will introduce a finite 
magnetization |14j . while in this work we are interested 
in equal population case. Hereafter we will only con- 
sider the polar type configuration. One can further show 
that the free energy is minimized when A sk and A t k 
have a relative phase 7r/2 which leads to A 2 2k = A k = 
(|A sk | 2 + |A tk | 2 )/2 and £ x , 2k = E k = y/£ + A k [TS]. 

Finally, we obtain the gap equation for A sk and A tk 
by minimizing the free energy functional Eq. Q as 

where £ = s, t, respectively. We note that in Eqs. |HJ 
A s and A t are not decoupled. In fact, they are coupled 
through £q which depends on both A s and At . The gap 
equation ^ can be solved numerically together with the 
number equation n = —dflo/dfi = j d 3 kn k /(27r) 3 where 
n k = 1 — tanh |^ and n = + n\. = kp/(3n 2 ) 
is the total density, with k-p the Fermi momentum of 
each spin component. Due to the azimuthal symme- 
try of the interaction potential, we shall consider the 
solution that A^ k only depends on k = |k| and cos6* k , 
and is independent of the azimuthal angle </> k of k. 
Furthermore, we can expand A^ k into different partial 
waves as A^ k = J2e A<^(/c)P£(cos#k) where Pe(x) is the 
Legrend polynomials and £ runs over even/odd integers 
for C = s/t. 

Results: In this problem, we can choose two dimen- 
sionless parameters D/a s and k-pD as control parameters 
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FIG. 2: (a) Ground state energy per particle as a function of 
fcpD at fixed D/ro marked in Fig. [ljb), D/a s < for dash 
green line with squares, D/a s = for blue solid line with 
circles and D/a s > for dash dotted red line with triangles, 
(b) Comparison of ground state energy between dipolar prob- 
lem and s-wave separable potential calculation (black crosses) 
with finite effective range. The inset: the comparison of or- 
der parameter, (c-d) The energy of pure singlet pairing (solid 
line), pure triplet pairing (dashed line) and mixed pairing 
(dash-dotted line), (c) is at the resonance point and (d) is at 
the BCS side. 

of the many-body physics. Since we focus on the regime 
nearby the first DIIR, there is one-to-one correspondence 
between D/a s and D/ro as shown in Fig. 1(b). And once 
D/ro is fixed, the other parameters of low-energy scat- 
tering process, such as s-wave effective range r e /D and 
p-wave scattering length a p /D are also fixed, as shown in 
Fig. 1(c) and (d). Another parameter k F D controls the 
ratio of dipolc length D, as well as the effective range r e , 
to the average inter-particle distance. The main results 
of this work are: 

(A) Nearby resonance where D/a s — 0, for small kpD, 
the system is in a pure singlet pairing phase. In this 
regime, despite that the dipolar interaction is anisotropic, 
the resulting order parameter is nearly isotropic. More- 
over, we show that in this regime the physics can be well 
captured by a simpler model, that is, a pure isotropic s- 
wave resonant interaction with a finite range correction. 

(B) As kpD increases, the system first undergoes a 
transition to a mixed parity phase with both non-zero sin- 
glet and triplet order parameters, which spontaneously 
breaks the time- reversal symmetry When kpD further 
increases, it experiences another phase transition to a 
pure triplet phase. A phase diagram across DIIR in terms 
of D/a s and kpD is constructed. 

Hereafter we shall explain these two results in details. 
(A) Small kpD Regime: We find that for small kpD 
the solution of gap equation Eqs. [8] gives pure singlet 



pairing, i.e. A s k ^ and A t k = 0. In Fig. [2]ja) we plot 
the ground state energy per particle E gs as a function 
of kpD for three different values of D/ro, as marked in 
Fig. 1(b). As noted above, once D/r is fixed, D/a s is 
also fixed, and the three curves in Fig. [2]ja) correspond 
to D/a s < 0, = and > 0, respectively. If we take the 
limit k F D — > 0, k F a s — > 0~ or 0+ for fl/a s < and 
> case, respectively. As shown in Fig. [2]ja), E gB ap- 
proaches free fermion case of 0.6-Ep in k F a s — > CP case, 
while E gs continuously decreases in k F a s — > + case be- 
cause all fermions form deep bound states. If D/a s = 0, 
as k-pD — > 0, E gs approaches 0.353EF, and very surpris- 
ingly, this coincides with the mean-field results of unitary 
Fermi gases with single channel s-wave contact potential. 
That is to say, in the regime of small k F D and nearby 
a DIIR, the pairing energy is on the same order of the 
Fermi energy as in a unitary Fermi gas at a magnetic 
Feshbach resonance, and thus, its transition temperature 
is expected to be as high as a unitary Fermi gas. 

Such behaviors are strongly in contrast to what one 
can obtain from Born approximation, where the sin- 
glet pairing is always exponentially small as kpD — > 0. 
Thus, these are unique features due to scattering res- 
onances. Furthermore, although the dipolar potential 
is anisotropic, the solution of order parameter is very 
isotropic in this regime. As shown in Fig. [3^ a), the order 
parameter is dominated by s-wave, and the contribution 
from d-w&ve and even higher waves are negligibly small. 
This is in fact consistent with the physical insight from 
two-body problem, where the scattering amplitude be- 
comes very isotropic at DIIR and the contributions of 
high partial waves become very small [U [TU] . 

Based on these observations, we can further show that 
the physics of this regime can be captured by a simpler 
model of an s-wave separable potential V£ k , = VoWkVOw, 
where Wk — 0(A— k) and a momentum cutoff A is intro- 
duced to implement a finite effective range. By solving 
two-body problem with VJ^., one finds a s — oo by fix- 
ing mVoA = — 27r 2 , which possess a finite effective range 
r e = 4/(ttA). Thus, we fix fc F /A = 1.157rfc F L>/4 so that 
the effective range from this s-wave separable model is 
fixed to the same value as the dipolar potential. As shown 
in Fig. [2^b), when kpD < 0.3, the ground state energy 
and the strength of order parameters obtained from this 
s-wave model agree very well with dipolar model at DIIR. 

The finite range correction can explain why energy 
increases as kpD increases in this regime. Because we 
can write an energy dependent scattering length a s (k) as 
l/a s (k) — l/a s — r e fc 2 /2, thus, with positive r e , a s (k) 
is negative at DIIR for finite k. At the Fermi energy, 
l/(fcpa s (fcF)) = l/(kpa s ) — kpr e /2. Thus, as kpD in- 
creases, kFr e increases and the attraction around the 
Fermi surface becomes weaker. Thus, the energy in- 
creases. Recently, several other papers have also studied 
finite range correction to an s-wave model |16) and our 
results are consistent with theirs. 
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FIG. 3: Partial wave distribution of superfluid order param- 
eter right at the first DIIR. (a)-(c) pure singlet pairing with 
fcpD = 0.1, 1.5 and 2.2 for (a), (b) and (c), respectively. 1 = 
for blue solid line; 1 = 2 for green dashed line; I = 4 for red 
dash-dotted line, (d) pure triplet pairing with kpD = 1.5: 
I =1, 3, 5 corresponding to the same symbols as t =0, 2, 4 in 
(a-c). 

Finally we show in Fig. [4] (b) and (c) that in the singlet 
phase, the system display standard behavior of BCS-BEC 
crossover, that is, the pairing gap continuously increases 
and the chemical potential continuously decreases as a 
function of l/(fcpa s ). 

(B) Large kpD Regime: As kpD increases, we first find 
that in the singlet pairing phase, the cZ-wave component 
and higher partial wave component gradually increase, 
as shown in Fig. [3|T>c). On the other hand, the en- 
ergy for pure triplet pairing (A t k ^ and A s k = 0) 
gradually decreases, this is because the triplet pairing 
is dominated by p-wave as shown in Fig. [3lc), and a p 
is proportional to D as shown in Fig. fTTcL) . There- 
fore, as kpD increases, the energy of the triplet pairing 
phase will become lower than the single pairing phase. 
At DIIR, we find that this energy level crossing occurs 
around fcj?D ~ 2. Moreover, in the regime when singlet 
and triplet paring energy are close, we find a new super- 
fluid state with nonzero A s k and A t k whose energy is 
lower than pure singlet and pure triplet paring, as shown 
in Fig. |2jc-d). Since A s and A t have opposite parity un- 
der the time-reversal transformation, thus, time-reversal 
symmetry will be spontaneously broken when they co- 
exist [T5] ■ As shown in Fig. [2jc-d), the parameter win- 
dow for mixed paring state decreases as one moves to the 
BCS side of the DIIR. Based on the energy comparison, 
a phase diagram is constructed as shown in Fig. |4^a) , in 
which two second order phase boundaries separate sin- 
glet, mixed and triplet paring phases. 

Hence, we have established results (A) and (B). 

Final Comments: At last, let us make two comments. 
Firstly, we notice that the Fock exchange energy will 
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s 

FIG. 4: (a) Schematic of the phase diagram for spin-| Fermi 
gases of polar molecules across a DIIR. (b-c), In the singlet 
pairing phase, A/£f (b) and u/£f (c) as a function of l/kpa a . 



deform the Fermi surface [T7H2D] and has influence on 
fermion pairing [3T] [52] . Such effect can be included by 
replacing the bare dispersion £k with ek + Sk where Ek is 
the Fock self-energy given as Ek = — J2q ^d(k ~ Q) n q/2- 
This equation for Ek needs to be solved together with 
Eqs. |8]self-consistently. We find that including this Fock 
exchange term does not lead to any qualitative change to 
the results discussed above. It only quantitatively shifts 
the phase boundary to lower value of kpD, which im- 
plies that the deformation of Fermi surface seems to favor 
triplet pairing over singlet paring. At small kpD < 0.5, 
even the quantitative correction to the ground state en- 
ergy from the Fock term is negligibly small. 

Secondly, we notice that there are several papers using 
Tk,k' = 47ra s /m + Vd(k — k') to study fermion pairing 
problem [15l [2TJ [22] , where they assume a s is originated 
from the short-range potential, and Vd is taken as the 
long-range dipolar potential. Such a treatment is also 
inappropriate near a DIIR. This is because although the 
interaction potential is a combination of the short-range 
and the long-range part, the scattering amplitude can not 
be taken as a sum of two contributions unless one takes 
Born approximation. Furthermore, in this treatment, if 
a s 3> D the contribution from dipolar part completely 
vanishes for all range of k-pD. Therefore, it can not cap- 
ture the finite range correction for small kpD and the 
phase transitions for large kpD. 

Acknowledgements. We thank Tin-Lun Ho, Georgy 
Shlyapnikov and Eugene Demler for helpful discussions. 
This work is supported by Tsinghua University Initia- 
tive Scientific Research Program. HZ is supported by 
NSFC under Grant No. 11004118 and No. 11174176, 
and NKBRSFC under Grant No. 2011CB921500. RQ is 
supported by NSFC under Grant No. 11104157. 

Note Added: During writing this manuscript, we be- 
come aware of another paper where the DIIR is modeled 
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in a completely different way [23. for studying many-body 
physics. The focus of these two papers are also different. 
Both the finite range effect and the phase transition be- 
tween singlet and triplet phases are not discussed there. 
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